arXiv:gr-qc/0512062v2 29 May 2006 


Modified Newtonian Dynamics and Induced gravity 
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Modified Newtonian dynamics, a successful alternative to the cosmic dark matter model, proposes 
that gravitational field deviates from the Newtonian law when the field strength g is weaker than 
a critical value go. We will show that the dynamics of MOND can be derived from an induced 
gravity model. New dynamics is shown to be compatible with the spatial deformation of scalar 
fields coupled to the system. Approximate solutions are shown explicitly for a simple toy model. 

PACS numbers: PACS numbers: 98.80.-k, 04.50.+h 


I. INTRODUCTION 

The rotation curve (RC) observations indicates that less than 10 % of the gravitational mass can be measured 
from the luminous part of spiral galaxies. This is the first evidence calling for the existence of un-known dark matter 
and dark energy. In the meantime, an alternative approach, Modified Newtonian dynamics (MOND) proposed by 
Milgrom JJ, has been shown to agree with many rotation curve observations QQ 

Milgrom argues that dark matter is redundant in the approach of MOND. The missing part was, instead, proposed 
to be derived from the conjecture that gravitational field deviates from the Newtonian 1/r 2 form when the field 
strength g is weaker than a critical value go ~ 0.9 x 10 -8 cm/s 2 Q. 

The phenomenological foundations for MOND are based on two observations: (1) flat asymptotic rotation curve 
(RC) is a common feature for many spiral galaxies, (2) the Tully-Fisher (TF) law, M ~ V a 3J is very successful 
explaining the relation between rotation velocity and luminosity in many spiral galaxies, a is shown to be close to 4. 
Note that first fact indicates that gravitational field g goes like 1/r asymptotically in the flat RC region. 

There have been effort trying to connect non-relativistic MOND theory with a relativistic versionfj. We will show 
that the dynamics of MOND could be due to the effect of scalar fields in an induced gravity model. It is known 
that scalar field can be used to induce mass for a massless fermion from a Yukawa coupling. The associated effective 
potential is also a source of dark energy. In fact, one of the original effect of the scalar field is to deform the definition 
of scale and distance in the Weyl invariant model. 

Indeed, Weyl US proposed that the invariant length scale should be defined as ds 2 = tyg^dx^dx 11 . It is easy 
to find that ds 2 is invariant under the local scale transformation <p(x) —> A 2 {x)(f>(x) and g^ —> A ~ 2 (x)g llv . The 
local length scale dr 2 = g^dx^dx” measured by any observers will then be deformed due to the scale chosen by the 
dimension two scalar field q 4. 

Therefore one of the effect of the scalar field is to deform the effective distance between two distant matters. One 
will show explicitly in this paper a possible resolution derived from a simple model with a spherically symmetric 
mass-distribution directed by similar effect. 

Einstein theory works with the latest CMB anisotropic survey in a manageable way. One would like to study if 
there is any way to induce the MOND under the Einstein-Hilbert framework with minimum modifications. Indeed, 
one is able to show that the dynamics of MOND could be due to the large-scale spatial inhomogeneity of the coupled 
scalar field. A specific example will be presented in this paper to show its effect in the large r limit where the effective 
potential V(r) goes like lnr in the absence of any dark matter. 

In section II, the idea of MOND will be briefly reviewed. The force law of MOND will also be rewritten in term of 
dimensionless field variables for convenience in this section. The effective potential for MOND will be shown in section 
III. In section IV, one shows that induced gravity model could provide a possible resolution to the force law of MOND 
with the help of two scalar fields ip and <f >. Two different approximate solutions to the field equations associated with 
a simple toy model are presented in section V. Some conclusions and remarks are drawn in section VI. 


II. MOND 

It was proposed that there exists a critical acceleration parameter go = 0.9 x 10 -8 cms -2 0 (Sanders 2002) charac¬ 
terizing the turning point of the effective power law associated with the gravitational field in MOND. Gravitational 
field of the following form was suggested 
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with a function /i( x ) considered as modified inertial. Here gjy is the Newtonian gravitational field produced by any 
sort of mass distribution. Milgrom suggests that the following function-form 
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agrees with the rotational curve of many spiral galaxies. 

It is interesting to find that there exists a length scale factor r o associated with our Milky Way via the following 
equation 


GM 0 
go = —— 
r o 


( 3 ) 


with Mq = 2 x 10 11 Hg roughly the order of total mass of the Milky Way. M s denotes the solar mass. As a result one 
can show that ro ~ 49700 ly which is roughly the size of the luminous part our Milky Way. 

Note that Mq is the typical mass of any galaxy with a size similar to our Milky Way. In fact, r o chosen here will 
only served as a convenient unit of length scale . Selecting different value of Mo and hence ro will not affect the 
physics presented in this paper. It turns out that writing all physical variables in dimensionless form will make it 
easier for us to focus on important physics. 

One will try to write the modified field strength g in a dimensionless form with the help of the parameters Mo, ro 
and go- The idea is to write M, r and g in units of M 0 , r 0 and go respecti vely. 

Milgrom argues that one can take the effective inertial g as g(x) = x/y /1 + x 2 . Therefore, one obtains 
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by dividing both sides of the equation pp by the critical parameter go- One can therefore write above equation as 


V 1 + g 2 
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with g and gw written in unit of go. This makes g and gw dimensionless from now on. 

Throughout this paper, we will focuss on the study of the system with a Newtonian attraction of the form g jy = 
Gm/r 2 for simplicity. This is the field strength at a radial distance r from a spherically distributed matter with total 
mass m. System with different mass distribution can be obtained straightforwardly. To summarize, one has chosen 
dimensionless scale according to following redefinitions: 
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For example, g' N = m'/r' 2 . One will suppress the superscript ' for convenience. Therefore, the modified field strength 
g becomes 


g 2 m 

s/TTg 2 ~ r 2 


( 7 ) 


One can further remove the parameter m by defining r c = roy/m and write r as a dimensionless coordinate variable 
in unit of r c instead of vq. 


r c y/mro 

As a result, one can write above equation in a very compact form 
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in terms of the dimensionless parameters g, gw, and r". Note that it is straightforward to restore all dimension 
parameters to evaluate any corresponding physical values. Suppressing the superscript " again for convenience, one 
can write it as 
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v/i+7 


( 10 ) 


Note that the dimensionless physical coordinate r" = r/r c = r/(y/m(r)ro) is in fact a function of r for a system 
with total mass distribution given by m(r). Gauss law implies that, for a spherically symmetric system, exterior mass 
throughout r > rq will not affect the gravitational field g(r < r\). From now on, one will assume m is a constant 
independent of the coordinate r for simplicity keeping in mind that the result is valid only for the exterior region of 
a spherical system. 

Once m(r) is not a constant, the generalization is still straightforward. One can simply take r" as co-moving 
coordinate. All physics can be re-derived by a proper coordinate transformation. 

Note that Eq. (11 011 can be solved directly to write g(r) as a function of r: 


9 = 
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It is apparent that g(r) goes like 1/r 2 at short distance scale where r<l. On the other hand, g{r) goes like 1/r 
at large distance scale where r>l Q. 


III. EFFECTIVE POTENTIAL OF MOND 


In order to take a close look at the changing pattern of g , one can split it into two different parts: Jj] 

2V2r 2 , y/l + Vl + 4r 4 

9 — 9i + 9s — -~ 7 - H- 7 = 2 rrr^rT' ( 12 ) 

\/l+4r 4 A/l + \/l + 4r^ V2r 2 Vl + 4r 4 

Note that the first term is gi that goes like 1/r when r > 1, while the second term is g s that goes like 1/r 2 when 
r < 1. Therefore it is easy to see that gi and g s represents the long distance and short distance field strength of the 
g respectively. One is hoping that successful separation of g may help shedding light to the search of the underlying 
theory. 

One can integrate g for the effective potential. After some algebra, one can show that the effective potential 

X T* pT* cT 

x gi ■ dr = J gidr and V s = J g s dr can be evaluated directly to give 


and 
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One can verify directly that V/ = gi and V/ = g s and prove that above equations are indeed correct up to an irrelevant 
integration constant. 

In fact, it is difficult to specify this constant of integration in the conventional approach which take V(r —> 00 ) —> 0. 
This is because the effective potential in fact diverges at spatial infinity due to the logarithm behavior of the dominating 
2D-like potential. 

One remarks that the potential V(r) derived here remain valid in the co-moving coordinate chosen as r/r c which 
depends on the mass content of the spherically symmetric system. 


IV. INDUCED GRAVITY THEORY AND MOND 

One is looking for a theory that goes to Newtonian theory in the small r region and reproduce the g ~ 1/r effect 
in the large r region. The dynamics proposed by MOND gives a specific force law shown earlier in section II. In 
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particular, the induced gravity (or equivalently the Brans-Dicke theory) provides a nice framework for our purpose. 
Note that induced gravity model, similar to the Brans-Dicke theory, proposes that the the gravitational constant is a 
dynamical variable given by the vacuum expectation value of a scalar field (p derivable from the Lagrangian Q 

^ = ( 15 ) 

Here R is the scalar curvature and <p is a scalar field producing a space-time dependent gravitational constant. In 
addition, W{cp) is a spontaneously symmetry broken effective potential coupled to the system. In addition to the 
Lagrangian of scalar field there is another matter Lagrangian 

C-m = C M (m,ip) ( 16 ) 

with a modified function ip(x) introduced here to induce the MOND effect. One will assume that the variation of 
C M with respect to the metric will contribute a term %’T^l in the field equation accounting for the matter effect. 
Here ip is treated as an auxiliary field without any dynamics. The main purpose of this auxiliary field is to introduce 
a distortion of Newtonian potential in the Newtonian limit. 

The field equations of the induced gravity theory are known to be (gj 

(pG „„ = -8 t ^T^-^D^D^-^D^cP) 

- {D^D^cp - g^D 2 p) + Aug^W ( 17 ) 

d ^ = 3T '-J T r + Sry a * w - 2W \ < 18 > 

with T^l the energy momentum tensor associated with the matter and ip another scalar field coupled to the mass 
term of the constituent matter. Note that the variation equation of <p will produce a term proportional to the scalar 
curvature R. This R proportional term can be eliminated by the R term obtained from taking the trace of the Eq. 
EJ. The final result is exactly the Eq. m shown above. 

Physical interpretation of ip field: 

Note that mass of a fermion is induced from a Yukawa coupling term, e.g. bip'b for a baryon b. The mass of the 
baryon b is thus induced by the vacuum expectation value of the coupled scalar field < ip' >—> mb- Assuming that 
rrib(x) =< ip'{x) > is spatial dependent, one can in principle induce a very different mass effect for the system in 
coherent to the spatial dependent scalar field (p{x) introduced earlier. To be more specific, the inhomogeneous nib(x) 
represents the local mass deformation in our approach. 

The overall effect can be integrated in order to produce global deformation on the Newtonian potential V{x,m). 
In our approach, one assumes that the deformation function ip(r) coupled to the matter energy-momentum tensor 
represents collectively the total effect of the inhomogeneity distortion of MOND potential V{x). In other words, ip is 
assumed to provide the collective deformation of the entire system accommodating the MOND potential as a physical 
resolution. In addition, one assumes that the mass generating scalar field ip' is an auxiliary field without a kinetic 
term. Moreover, one would like to impose the equation 

ipV 0 = <pV (19) 

as the auxiliary constraint for ip. Note that what one imposes here is a relation between ip and V instead of relating 
ip and (p. Indeed, one can see from the structure of the field equation GO) that only the ratio ip/(p coupled to 
matters. Indeed the field equation reads G^ = [ ip/(p]T+ ■ ■ ■ with the term [ip/(p]T^l serves as the generalized 
energy momentum tensor. Therefore, one can parameterize ip differently, but the final result will remain the same. 
By all means, one finds that the constraint (TTTT1) is the best way to relate ip and V. This constraint introduced here 
is similar to the effect of modified inertial introduced by the scalar-tensor theory 0 ]. By all means, it requires a 
constraint to be imposed by hand in order to reproduce the dynamics of MOND in an exact form. Hopefully, the 
approach shown here will shed light to the finding of a more realistic approach. 

In fact, the constraint one introduced here is equivalent to setting ip/cp = V/Vq such that the system will become 
Newtonian in the limit V —> Vq. To be more specific, the constraint is in fact introduced via the ration of (p/ip . Bring 

ip/(p to the right hand side of above equations, one can write G = [ip/4>\T^• ■ ■ as (p'G^ = [cp / ip]G^ = T^-\ -. 

This is equivalent to the induced gravity theory without a ip field. One can introduce the constraint via a re-scaling 
of the scalar field (p’ = <p/ip in the same induced gravity theory. In this approach, one can as well say that the 
constraint is imposed on the choice of scalar field ip' = (p/ip by choosing a proper gauge of scale transformation |f|. 
And hopefully, the constraint cp' = Vq/V (equivalent to the constraint 1191) 1 governing the deformation information of 
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the Newtonian force law will also lead one from Newtonian theory to the dynamics proposed by MOND. In short, one 
can either view the constraint m as the collective mass function ip or the deformation of the Newtonian constant G 
prescribed by <p>'. 

Note however that the field equations for these two different approaches will be slightly different sue to the kinetic 
term of the scaled phi field. In a moment we will try to solve this model in the limit where kinetic term of <p is omitted. 
The field equations of the scalar field will be the same in this limit. Therefore, it does not matter which viewpoints 
one decide to take in the toy model we will study later. The view of the collective mass deformation is however a 
quite promising idea, therefore, one will stick to Eq.s (tTTITHl) throughout this paper. In short, the auxiliary field ip 
introduced here is not a dynamical field. In addition, the effect of ip has already been absorbed into the collective 
mass term included in ipT^ u . 

The classical Newtonian field equation can be obtained from the time-time component of the Einstein equation in 
the Newtonian limit with the metric identification 


goo = — 1 — 2V. 


( 20 ) 


Indeed, the geodesic equation 


becomes 


d 2 x M M dx a dx 13 
d S 2 ais ds ds 



( 21 ) 


( 22 ) 


in the Newtonian limit. Note that one can write V = tp(r)V o for convenience with V 2 Vo = 47 t p the Newtonian 
potential for a spherically symmetric system. Therefore ip = V/Vq represents the deformation factor signifying the 
deviation of the gravitational potential V as compared to the Newtonian potential Vo- One would like to study if 
there exists a consistent solution in an induced gravity model that will accommodate the gravitational potential of 
the form V(r) = (V) + V),) given by Eq. 11311411 . 

Note that the time-time component of the Einstein equation can be shown to be 

+ V)[<p" + ^V] - [Vo(ip" + ^V) + 2y 0 V]</> - (p^y 2 Vo 
+4nipp = -4w(i + E)W-w(i + v )tpj (23) 


with uj a coupling constant. 

Since one has V 2 Vo = 47 rp, the constraint ipV 0 = <pV implying ip = (pip can be used to write Eq. iPl) as 


(- + V)[(p" + -<p'] - [Vo{ip" + -ip') + 2 Vop'](p 

= -4w(± + V)W-u(± + V)^ (24) 

In addition, one has 

r+(- r +vyu = - m ( 25 ) 

from Eq. 11711811 . 

Therefore one has a set of equation with two correlated ODEs and an unknown scalar potential W(<p) to be dealt 
with. In principle, one can choose some appropriate, possibly exotic, combination of scalar potential W such that the 
equations (23-24) do accommodate consistent solution for (p. One should first solve the second Eq. (24) to find the 
solution of <p{r). Substituting this solution back to Eq. (23), one should be able to find ip(r). 

In practice, what one actually is doing is: (i) given any prescribed form of Newtonian potential V (r) and insert this 
function to equations 1241251) . (ii) one can in principle solve for a consistent set of solution (</>(r), W(r)). Finally, one 
can reconstruct W(r) as a function of < p(r). This in principle will give us a solution for W(<p). One may however end 
up with a very complicate expression for the scalar potential W. 

To be more specifically, one can write the scalar potential as 

-(- + V)[(p" + ~(p'] + [Vo (<p" + -ip') + 2Vq(p'](P 
2 r r 


M- + v ) 


w = 


(26) 
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from Eq. In addition, one can put Eq. lf25l as 


v'W 


87T 


3 -|- 2 uj 


Wp- 


8 oj 


3 T 


\4>w' - 2 <f>'W] 


(27) 


by multiplying it with <j>'. Eliminating W' and W, with the help of Eq. CB, one can write Eq. as a differential 
equation of </> independent of W with a given V(x) and p(x). This equation can be solved to give a formal expression 
of <j>(r). Insert this solution of </>(r) back to Eq. 12611 . one can hence write the scalar potential W(r) as a function of r. 
Since we know the form of (j>(r), one can invert the function to find r = r((f>). Therefore, one should be able to write 
W((f>(r)) as a functional of <f>(r) straightforwardly. One will show later as an example that the process does work in 
the large r limit where the gravitational force become dramatically deformed to the MOND limit g ~ 1/r. 

For our purpose, one needs to know if the set of field equations accommodates arbitrarily specified (p(r) as a 
consistent solution with a properly chosen W. One has two variables (f> and W{<j>) at our disposal. Therefore, the 
answer turns out to be yes. Indeed, one can always solve <j>(ip) as a function of <p with any given <p(r) and W (</>) from 
Eq. (23). The additional Eq. (24) would then require a very special form of W in order to have the needed function 
</(</?(?')) as a consistent solution. 

Since consistent solution with any given ip(r) can be in principle made possible with certain properly chosen potential 
coupled to the system. Given a carefully chosen potential W, one is naturally lead to the desired deformation factor 
<p(r) that induce MOND as an alternative theory via the identity V(r) = <p(r)Vo(r). 

Note that one assumes that the scalar potential W could be of very complicated origin. It might have to do 
with some complicated gravitational field interaction similar to the temperature-dependent effective potential. One 
is unable to offer a solution here. Hopefully the approach shown here may shed a light to a better understanding of 
the physical origin of MOND. 

For a simple demonstration, one will show explicitly in next section how to obtain a desired asymptotic solution 
with a toy model. 


V. A TOY MODEL 

One will try to show that the system of equations does accommodate consistent set of solution for cf> and 

W ( <f >) in the large r region in this section with the help of a simple toy model. For simplicity, one will consider 
the limit u —> 0. In order to keep the effective potential term W\ in our toy model, one also has to assume that 
W\ =W/uj. Equivalently, one is trying to remove the effect of the kinetic term from the field equations such that the 
scalar field 4> becomes an auxiliary field without dynamical term. To be more specific, one is considering the scalar 
Lagrangian given by 


C^-^R-WW. (28) 

We will try to find models which enable a MOND solution in large r region where no matter is present. In this case, 
the Newtonian potential Vo —> 1/r according to the Gauss law. 

In constructing the model, what one actually did is that: (i) let V —► ln(r/ri) be an approximate solution in this 
limit, (ii) insert this solution to equation mm and solve for a consistent set of solution ((j>(r),W(r)). Finally, one 
can reconstruct W(r) as a function of <j>(r). 


A. weakly coupled V 


In the region V <C 1, W <C <j>" and (j)W' <C 4>", the field equations become: 

0" + - 4 ? - 24V '4 

r 

4" + (- + v'W - o 

r 


(29) 

(30) 


with a potential W is either 0 or close nothing in this region. Note that one has used the identification V = ipV o and 
Vo = —1/r in this limit. 
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It is straightforward to show that the following solutions 


0 - 

00 
r 2 

(31) 

V - 

4 In — 

(32) 


solve Eq. (I29I30H in a consistent way. Here the parameter n represents a local re-scaling of the zero point of the 
gravitational potential. The constant — In n will not affect the force law Eq. El- 

In addition, the constraint V l implies that this solution is good for r close to r\. For r\ ~ ro, the solution does 
reproduce the 1 /r force law beyond the luminous region of spiral galaxies. It was shown in Ref. *7j that V —> lnr 
quickly when r > r c . Hence region close to r c is a perfect domain for our assumption. 

Note that the assumption W <C ft' and ftW' ft' implies that the scalar potential W{ft) = IE(0o/r 2 ) ft' ~ 1/r 4 
and W'(ft) = IE'(0o/r 2 ) <C 1/r 4 in this region. For example, W = A0 4 is a perfect candidate for the scalar potential 
near the region r ~ r±. If ri is bigger than the galactic size, this solution provides a reasonable and good prediction 
for many RC observations. 

One knows that the scalar potential W may be affected by the detailed dynamics of the system. For example, it 
will depend on the temperature of the system |fj. Therefore, it is also natural to expect a potential correlated to the 
scale length r in a large-scale astrophysical system. Once the effective scalar potential becomes negligible in a domain 
near and beyond our galactic size, the toy model discussed here may provide a good resolution for the new dynamics 
in that region. 


B. large r region 

One will show that if the scalar potential is given by 

W --^r exp[—2</>/ ft] (33) 

in the large r region, one will be able to reproduce a gravitational potential V —> lnr asymptotically in the large r 
region. 

Note that we are assuming that the weak field approximation goo —> — 1 — 2V and the resulting field equations 
(23-24) still hold in the large r limit even the assumption V <C 1 is not valid in the large r limit. 

There are two main reasons for this approach. Firstly, MOND approach also have similar problem with the 
asymptotic divergent potential V that goes like In r in the large r region. One has to deal with the potential difficult 
region in any case. If the induced gravity model considered here has anything to with the real physics, there must, 
hopefully, exist certain sort of normalization process, similar to the re-normalization theory, holding Eq. mm 
valid in the asymptotic region. Secondly, the purpose dealing with this toy model is simply to demonstrate how it is 
possible to derive the desired large scale (or weak field) gravitational potential V consistently, at least in the large r 
region. 

The field equations become 


Vft' + -Vft ~ [V 0 (<p" + -ft) + 2Voft](j) - AVW (34) 

r v 

ft' + (- + V')ft ~ - 2W] 

V o 

~ 200 exp[—20/0 o ] - 8W. (35) 

It is easy to show that Eq. does have an asymptotic approximate-solution 0 —> 0o In r by ignoring the negligible 
terms proportional to W. One can readily show that the W term is indeed negligible when 0 —> 0o In r in Eq. on. 
Substituting this asymptotic solution 0 into Eq. (El. one can derive the following equation 

lnr ft' ~ -^tp (36) 

Note that one has also ignored the negligible IE-dependent term in Eq. El- One can hence easily show that the 
solution to above equation is 


ip — r lnr 


(37) 






which implies the asymptotic solution V —> In r in the large r region. 

Therefore, one shows that it is indeed possible to derive the asymptotic solution with a toy model given by the 
action G3). In principle, it is possible to reproduce any physical deformation of the gravitational potential one desires 
with the help of some properly chosen scalar potential. 

Note that one in fact expects that V —> In r and hence ip —> —rlnr as give by Eq. Therefore, one inserts 

this solution back to the held equation and try to find potential W which works along with the held equations. As a 
result, one hnds that W ~ —3<^§ exp[—2</i)/<^>o]/8<?i> is the effective potential one is looking for. 


VI. CONCLUSION 

The force law of MOND is reviewed and rewritten in term of dimensionless held variables for a spherically symmetric 
system in this paper. 

One of the effect of the scalar held is to deform the effective distance between two distant matters. The scalar held 
ip coupled to the matter energy momentum tensor is proposed to represent the spatial inhomogeneous effect of the 
distance deformation in a collective way. One has also shown explicitly in this paper a possible resolution derived 
from a simple model incorporated with a spherically symmetric mass-distribution generated by similar effect. 

Indeed, one shows that the dynamics of MOND could be due to the large-scale spatial inhomogeneity of the coupled 
scalar held. A specihc example with a toy model is presented in this paper to demonstrate its effect in two different 
regions: (i) one presents an approximate solution in the limit V ~ ln(r/ri) -C 1 which holds in the region near r = r\. 
(ii) approximate solution is obtained in the large r limit. Note that the effective potential V(r) goes like lnr in the 
absence of any dark matter in this toy model. Indeed, 

Eq.s 1531251 are in principle solvable with some properly chosen potential even one is unable to obtain an analytical 
solution at the moment. The method shown in this paper also applies to different force law, with a different <p{r), 
that may properly describe the RC with or without any sort of dark matter. The constraint (pi) is introduced as an 
auxiliary constraint. One does not know the physical origin of this constraint other than the consistency in dimension. 
It deserves more attention studying the physical origin of this constraint. 
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